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Electron transmission through nanosystems is blocked if there are no states connecting the left 
and the right reservoir. Electron-electron scattering can lift this blockade and we show that this 
feature can be conveniently implemented by considering a transport model based on many-particle 
states. We discuss typical signatures of this phenomena, such as the presence of a current signal for 
a finite bias window. 
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I. INTRODUCTION 



II. THE MODEL 



A. The System 



Quantum dots constitute an excellent testbed for 
transport through general nanosystems, where the local 
density of states is dominated by discrete localized levels. 
The key points are conduction quantization^ due to the 
discreteness of levels, Coulomb Blockade due to electron 
repulsion^ and the interplay between resonant tunnel- 
ing and charging in double dot structures^ In this work 
we consider a further issue, the transport by electron- 
electron scattering. 



Electron-electron scattering is not included in stan- 
dard transmission models^ where the Coulomb inter- 
action is taken into account by a mean-field approach 
frequently including exchange-correlation interactions as 
well. Within such models electron transport strongly de- 
pends on the presence of states in the system connecting 
both leads £ Here we show that electron-electron scatter- 
ing allows for additional transport channels and that it 
can be consistently implemented using a many-particle 
basis following the concepts developed in Refs. 001111, 

nam 



While in double-dot structures, each dot has direct ac- 
cess to a reservoir with a continuous level density, the 
situation is essentially different in triple-dot structures,— 
where the states in the central dot only couple to discrete 
states in the neighboring dots. Thus the properties of 
these states are far more sensitive to scattering events, 
which may essentially determine the transport through 
the structure. This is precisely the situation depicted in 
Fig. [TJ Here the upper level 4 of the middle dot can be 
filled from the left lead by resonant tunneling via level 1, 
while its lower level can be emptied into the right lead 
by resonant tunneling via level 6. Thus the current is 
very sensitive to scattering between level 4 and level 3. 
In this work we restrict to electron-electron scattering, 
which is appropriate if the phonon energies do not match 
the transition energy. 



The system depicted in Fig. [TJ is described by the 
Hamiltonian H = Fdots + -ffieads: 



H, 



dots 



a}a,i 



(1) 



refers to the dot region where a\ {a{) is the creation (an- 
nihilation) operator for the i'th state. Assuming that 
the states in the individual dots are strongly localized, 
only states in dots next to each other couple and we re- 
strict to those couplings fiy depicted in Fig. [TJ For the 
Coulomb part H ee we neglect interactions between the 
leads and the dots as well as interactions between next- 
nearest neighboring dots. Then we obtain 

H cc = U(a\a\a,2a,i + a^a^a^a^, + a^aQaQa^) 

+ U n {a\a\ + a\a 2 + a\a 5 + ala e )(ala 3 + a\an) 

+ {U sc alalaia 4 + U sc alala 3 a e + h.c.} . (2) 

Here U and U n are the matrix elements of the stan- 
dard Coulomb repulsion between states located in the 
same and neighboring dots, respectively. U sc describes 
Coulomb scattering between different states, which is the 
central issue of this work.— 

Finally, the Hamiltonian of the leads and their coupling 
to the dots reads: 



H 



leads 



^2 Euc kl cu 



ikl 



tiua\cu + tl u cl,ai (3) 




FIG. 1: The triple-dot system considered. The bold lines 
represent the simplified system. 
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Ei = 40 


Ei = 60 


£3 = 20 


£4 = 40 


E 5 =0 E 6 = 20 


f2l4 = ^36 = ^23 = f^45 = 0.1 


«13 = ^35 = 0.05 


r L = r R = o.i 


/Lti = 50 


HR = 10 


SI24 = ^46 = —0.2 


U ac = -0.2 


U = 10 


t/n =3 


fc s T = 2 


W = 400 



TABLE I: Parameters used if not stated otherwise. They 
refer to the modulated nanowire discussed in Sec. Ill Bl but 
have been rounded off for an easier recognition of scales in the 
plots. All energies are in meV. In the simplified system we 
neglect level 5 and set f2i3 = Q24 = ^23 = — U = U n — 0. 



Here £ = L,R denotes the left and right lead, respec- 
tively. The energies Em in lead t provide a continuum of 
states (labeled by k). We assume that the correspond- 
ing density of states has the constant value ge in the 
energy range —W<E<W and is zero otherwise. 
Disregarding the fc-dependence of the tunneling matrix 
elements t ike we set t lkL = t 2 kh = ^T L /(2Trg L ) and 

t5kR = t^kR = \/Tr/ (^9r)- These are the transitions 
sketched in Fig. [1] All other tunneling matrix elements 
are neglected in Eq. ([3]). Throughout this work we re- 
strict to a single spin direction for simplicity. 



B. Parameters 

For specific calculations we use the parameters of Ta- 
ble[T]unless stated otherwise. They relate to an InAs/InP 
modulated nanowire structure similar to the structures of 
Refs. H^flol . We assume three In As wells with a thickness 
of 40 nm, which are separated by 3 nm thick InP barriers. 
The outer barriers are assumed to be 1.5 nm thick. As 
we are only interested in order of magnitude estimates 
we choose the simple one-band envelope function model, 
with Dirichlct boundary condition on the outside of the 
wire. In addition, we assume that the wire is cylindrical 
with radius R = 20 nm which enables us to reduce the 
problem to a one dimensional problem by using cylin- 
drical coordinates, i.e., the single particle Hamiltonian is 
given by 



H=\- 



d h 2 d 



dz 2m c s dz 2m c ffi? 2 



Jin 



+ V c (z)] (4) 



where ji n is the n'th zero of the Besscl function J;. m e g 
is the effective mass function, V e g is the conduction band 
edge function (they are stepwise constant). In the pre- 
vious section we assumed that states in individual dots 
are strongly localized. One way of achieving that is to 
use Wannier states ^i(z) for individual dots, assuming a 
periodic repetition of the structure. Using the masses 
™inAs = 0.026m e and mi n p = 0.08m e , where m e is 
the free electron mass, and a conduction band offset of 
O.QeV— we get an energy difference between the ground 
state and the first excited state of E% — E\ = 21 meV 
for a given I and n. The excitation energy for the ra- 
dial modes 2 m it R: i Oil — Joi) = ^9 meV is larger and 



thus these radial modes can be neglected. (In addition, 
the coupling between states of different radial symmetry 
should be small.) The couplings f^- are evaluated fol- 
lowing Sec. 2.3 of Ref. for a bias drop of 20 meV per 
period. 

For the coupling to the leads we use the estimate^ 



2dT? 



iL/R 



h^2{E t - E r )/m ln 



(5) 



As 



where is the coupling element between Wannier state 
i in neighboring dots for a barrier width of 1.5 nm (the 
outer barrier) and E r is the sum of the conduction band 
edge of InAs and the radial confinement energy. 
In general the Coulomb interaction is described by 
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Hee U % jkia\a^a k ai 

ijkl 



with 

Uijki = 



d 3 r 



dV 



e 2 ^rV,(rM(r> fc (r') 
47re r e |r - r'| 



(0) 



(7) 



Commonly, one focuses on the direct interaction of two 
states, where i = I and j = k. Taking into account the 
normalization of the wave functions, we can estimate 



Ui 



47re r eorf 



(8) 



where d is the average distance between the particle den- 
sities. Using d « 10 nm, if the states i and j are within 
the same dot, and d m 40 nm, if the states i and j are 
within adjacent dots, we obtain the values for U and U n 
given in Table HI respectively, for e r w 13. 

The key scattering element U sc corresponds to U^u- 
As the states ^(r) and 954 (r) [as well as ^(r) and y>i(r)] 
are orthogonal, one cannot approximate the l/|r — r'| 
potential by a constant value as in the case of the direct 
interaction discussed above. Instead a dipole expansion 
is possible providing 



U: 



2^21^34 



3214 



47re r eo d 3 



(9) 



where d as 43 nm is the distance between the centers 
of neighboring quantum dots. The z-matrix elements 
Zij = J dztp*(z)zifj(z) are evaluated for the Wannier 
functions, providing Z12 = Z34 = —8 nm, which gives the 
value in Table |U As U sc -C U,U n , it is usually neglected. 
However, here we show that it can have an crucial impact 
on the transport. 



C. Transport approach 

For our calculations we use a basis of many-particle 
states j a), \b), . . ., which diagonalize the dot Hamiltonian 
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-ffdots includin g th e Coulomb interaction. Using the ap- 
proach of Ref. lid , but only including first-order transi- 
tion processes between the leads and the dot region, the 
following rate equations (first-order von Neumann ap- 
proach, see also Ref. [HI) can be derived for the reduced 
density matrix of the dot Ww = Tr{(6|/5|6')} where the 
trace is taken over all lead states {k£}: 



\h—w bb > = (E b - E b i)w w 

+ Yfo»W<fc a M - Aa(ke)T b *, a (k£)] 

+ ^T[T* b (k£)cl) c b>(k£) - fc b {M)T cV {kl)} , 

with T ba (kl) = Y,i Uu{b\a\ \a) and 
T cb ,(M)f e (E k ) 



k D T=0.5 meV 



b cb {k£) =J2 



Eh — E r 



iO" 



-w b > b 



E 1 c' 



T c , b {kl)[l- f t {E k )] 
E k -E c + E b + 10+ 



(11) 



Here f e (E) = (1 + e (E-^)/k b Tyi ig the Fermi digtri _ 
bution for the lead £ with electrochemical potential pig. 
The current from lead £ into the sample is given by 
h = J2cbJe( cb )i where 



Jt{cb) 



(12) 



is the part of the current associated with transitions be- 
tween states b and c within the dot i 10 ' 18 We disregard 
the sign of the electron charge e, so that the sign of the 
electrical current equals the sign of the particle current. 

It should be noted that we obtain the Pauli Master 
Equation^ if we neglect the off-diagonal elements of the 
density matrix. However, this approximation is only rea- 
sonable as long as the spacing between the many-particle 
energies is large compared to the contact couplings T^S- 
This is not the case for the systems considered here, so 
the full set of equations is needed. 



III. RESULTS FOR THE SIMPLIFIED SYSTEM 

At first we study a simplified system where we neglect 
state 5 and all interdot tunneling processes except for 
£1,23 and £736 which are in resonance. This corresponds 
to the thick lines and arrows in Fig. [T] In order to avoid 
complications due to Coulomb charging, we set U n = 
U = 0, thus focusing on the scattering via U sc . 

In Fig. [5] we show the current as a function of the 
left Fermi level /jl. There is no current until the left 
Fermi level comes within the vicinity of the ground-state 
of the first dot (£i = 40 meV). At this point electrons 
start to flow from the left lead into this state and further 
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FIG. 2: Currents as a function of fiL and £2 (subfigure) 
Other parameters from Table |T] (simplified system) . 



into the excited state of the second dot. If both states 
1 and 4 are occupied, the Coulomb scattering via f7 sc is 
possible. This process transfers one electron from level 
1 to level 2 and a second electron from level 4 to level 
3, which can subsequently reach the right lead via level 
6. Thus Coulomb scattering establishes a transport path 
through the nanosystem. However, when the left Fermi 
level comes into the vicinity of the excited state of the 
first dot (E2 = 60 meV) electrons will start to occupy 
this state. This causes a decrease in the current (see 
Fig. \2§ as Pauli blocking hinders the scattering process 
addressed above. Likewise the temperature dependence 
essentially follows the probability 

F = /i(£i)/ £ (£4)[l - f L (E 2 m - f R (Es)] (13) 

to find states 1 and 4 occupied while states 2 and 3 are 
empty. The relevance of level E2 for the transport is 
further demonstrated in the subfigure of Fig. [21 showing 
that current only flows through the triple-dot structure 
if E2 — Ei w £4 — £3, where the Coulomb scattering is 
energetically allowed. 

This presence of current enhancement in a finite bias 
window AfiL matching the energy transfer A£ by the 
scattering process is the characteristic signal of electron 
transport by Coulomb scattering. This Pauli-blocking 
of the scattering from level 4 to level 3 by occupation 
of the further level 2 does not appear for other inelastic 
scattering mechanisms such as phonon scattering. 



IV. DESCRIPTION BY SCATTERING 

The description based on scattering given above be- 
comes quantitative if the Coulomb scattering is the lim- 
iting process for transport through the device, i.e., if U sc 
is significantly smaller than O14 and r^g. For this reason, 
we have performed calculations for the increased values 
f2i4 = il36 = 1 meV, see Fig. [3] The strong coupling 
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FIG. 3: Current as a function of U sc with fli4 = i}-^ = 1 meV 
and E2 = 59 meV. The two subfigures show the current for 
U sc = 0.01 meV as a function of \xl and E2, respectively. 
Other parameters from Table U (simplified system). The solid 
lines are calculated by the first-order von Neumann approach 
and the dots depict Jl — eRti, where R tr is the transition 
rate, Eq. (|14p . for the electron-electron scattering process. 



between the states 3 and 6 yields a bonding and an anti- 
bonding state, with energies 20 =F 1 meV. Therefore the 
resonance condition for Coulomb scattering is now satis- 
fied at E2 = 59 meV and E2 — 61 meV (not shown) as 
displayed in the right subfigure of Fig. [3] 

Fermi's golden rule provides us with the transition rate 
by Coulomb scattering into the anti-bonding state be- 
tween 3 and 6: 



■Rtr 



U 2 



eff 



2 lA F ( 14 ) 



Here we have replaced the energy-conserving ^-function 
by a Lorentzian, representing life-time broadening due to 
the coupling to leads. r or T = 2r^ + is the sum of 

broadenings for the individual states: Tl for the levels 
1 and 2, and Tr/2 for the anti-bonding combination of 
3 and 6. Fig. [3] shows that Fermi's golden rule provides 
a full quantitative description for small U sc . However, 
for larger values of U sc , this simple reasoning, based on 
single-particle states, fails. In particular, the width of 
the current peak becomes much broader than the simple 
life-time broadening T e g (see subfigure of Fig. , which 
makes it easier to observe the effect in a real system with 
imprecise control over the level energies. 



V. DESCRIPTION BY MANY-PARTICLE 
STATES 

Now we want to sketch, how this scattering-induced 
transport emerges within a basis of many-particle states, 
which takes into account the entire Coulomb interaction. 
For the parameters of Table [J the anti-symmctrized two- 
particle product states |1, 4), |2, 3), and |2, 6) all have the 
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FIG. 4: a) Energies of selected two-particle states as a func- 
tion of E2- b) Coefficients of state |c). Other parameters from 
Table [T] (simplified system) . 



J L (cb)=1.27nA 



Jlb> 



□ 

;J R (c'b)=-0.73 nA 




:J L (cb')=-1.46nA 



J L (c'b')=0.73 nA 



/ 



Ilc'> 
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□ 



Jla> □ J R (b'a)=-0.73 nA 



FIG. 5: Diagrammatic representation of a circle of transport 
processes involving the many-particle states \b), \c), \b'), \a), 
and jc). The two sets of double lines are the bonding and anti- 
bonding combinations of the states 1, 4 and 3, 6, respectively. 



same sum of single-particle energies Ei + Ej = 80 meV. 
They couple to each other due to the matrix elements 
U sc and ^36, resulting in the three many-particle states 
depicted in Fig. [5] For E 2 ~ 60 meV, the three states 
are highly entangled and we focus in the following on 
one of these entangled states, denoted by |c). This state 
contributes to a circle of transitions between different 
many-particle states depicted in Fig. [5) The state |c) can 
be reached by tunneling of an electron from the left lead 
into the state \b) (process from upper left to upper right). 
Here \b) is the binding one-particle state combining levels 
1 and 4. By removing an electron towards the left lead (at 
a higher energy than before), the state |c) decays to the 
one-particle state the binding state combining levels 
3 and 6. Then the original state \b) is restored by one 
electron tunneling from the state \b') to the right lead 
and one electron tunneling into the state \b) from the 
left lead, which can happen in two different sequential 
orders. The key issue for the existence of this circle is 
the presence of the entangled state |c), which enables 
the transition between \b) and \b') via two single-electron 
tunneling processes. Therefore the current drops, if the 
product states |1,4) and (|2,3) + |2, 6>)/v / 2 arc detuned 
by varying £2- 

The currents JL/n(cb) in Fig. O denote the contribu- 
tion of the transitions b <-> c between the corresponding 
many-particle states to the current from the left/right 
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FIG. 6: Current as a function of (II for different temperatures. 
Other parameters from Table [I] In the subfigure the single- 
particle energies E\ = 43 meV, E2 = 63 meV, E3 = 20 meV, 
Ei — 40 meV, Es = —4 meV, and E% = 16 meV are used. 



observe a current signal in a finite region of (1l , which is 
a key feature of the current induced by electron-electron 
scattering (the current is below 0.006 nA if U sc = is 
used). However, contrary to the simplified system in 
Fig. H the peak current increases with temperature for 
the full system and is much weaker. This is due to the 
presence of an electron in state 5 which breaks the align- 
ment between the levels 3 and 6 by Coulomb repulsion. 
With increasing temperature, the probability for state 5 
to be empty increases and so does the current. In the 
subfigure of Fig. [6] we show results for a case where the 
single-particle energies have been modified to compen- 
sate for charging effects, which provides results similar 
to Fig. [21 In both cases we observe enhanced current in 
(multiple) finite bias windows Afi^ matching the energy 
transfer AE, which are however smeared out by temper- 
ature. This shows that the essential features of transport 
by electron-electron scattering are robust with respect to 
other electron-electron interaction mechanisms. 



lead into the system, respectively, as given in Eq. (|12[) . 
The magnitude of these currents corresponds to the tran- 
sition rate between the states. Fig. [5] shows that the in- 
going and outgoing rates partially balance for all states 
depicted. Nevertheless, there arc plenty of further tran- 
sitions, which make the full picture far more involved. In 
total this circle provides Jj, = 1.71 nA and Jr = —1.46 
nA, which constitutes only a part of the total current 
Jl = — Jr = 5.9 nA. The remaining part is carried by 
similar circles involving the other many-particle states 
as well as more complicated transitions which cannot be 
separated into circles that easily. 

Finally, note that the electrons enter the structure with 
E f=a Ei from the left contact and leave the structure with 
E k, E5 to the right contact as well as with E f=a £2 to 
the left contact. Thus there is no single transmission 
channel at a given energy as typical for the frequently 
used transmission models. 



VI. RESULTS FOR THE FULL SYSTEM 

In Fig. [5] we present results for the full system, i.e., 
using all parameters given in Table [J Like in Fig. [5] we 



VII. CONCLUSION 

We have shown that Coulomb scattering provides a 
current channel for transport through a triple-dot sys- 
tem. The mechanism holds for general nanosystems ex- 
hibiting two pairs of states with a similar level spacing 
AE. An example is the conduction through a macro 
molecule, where an appropriate chemical group appears 
twice. A typical signature is a current signal for a fi- 
nite bias window, matching the energy transfer AE. If 
the Coulomb scattering is the slowest transfer process in- 
volved, a simple description based on Fermi's golden rule 
is valid. Otherwise a systematic implementation is possi- 
ble within a basis of many-particle states, which reflects 
the total Coulomb interaction for the nanosystem. 
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